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An exact theory for the density of a one-dimensional Bose-Einstein condensate with hard core 
. particle interactions is developed in second quantization and applied to the scattering of the con- 

■ densate by a spatially periodic impulse potential. The boson problem is mapped onto a system of 
Q\ | free fermions obeying the Pauli exclusion principle to facilitate the calculation. The density exhibits 
t-H ■ a spatial focusing of the probability density as well as a periodic self-imaging in time, or Talbot 

» I ' effect. Furthermore, the transition from single particle to many body effects can be measured by 

observing the decay of the modulated condensate density pattern in time. The connection of these 

■ results to classical and atom optical phase gratings is made explicit. 
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I. INTRODUCTION 

The recent interest in neutral atom Bose-Einstein con- 
densates (T]J^] and in periodic effects ||-|^] in atom optics 
has stimulated research in those fields. Experiments are 
at the point now where the condensate can be manipu- 



q , lated U and used for precision measurement and atom 
optical studies. From a theoretical perspective much can 
be learned by applying the methods developed over the 
past four decades in condensed matter research to these 
Oh| fields. Similarly, new effects from a condensed matter 
viewpoint can be introduced by adopting the knowledge 
gained over the years in mathematical physics, optics, 
J> , and atomic physics. With this motivation we detail the 
theory of a one-dimensional condensate with hard-core 
s ! \ particle-particle interactions subjected to a spatially pe- 
riodic, impulsive interaction. The subsequent evolution 
of the system probes the onset of many-body effects rep- 
resented by the interactions. 

The results are summarized as follows. The conden- 
sate of interacting bosons can be mapped onto a system 
q ' of free fermions obeying the Pauli exclusion principle. 
The boson-boson interaction in one-dimension acts like 
the Fermi repulsion to create a parabolic band of occu- 
pied states. If we are interested in observables, such as 
O |. the density, which probe the diagonal density matrix el- 
ements of the system, the difference between the fermion 
• i-h , and boson commutation relations is irrelevant to our re- 
^ ■ suits. 

Using this fermion model, we show that an exact ex- 
pression for the density after the periodic interaction can 
be derived for the case of a pulse with duration r such 
that max[efc F r, e q r] <C h, where the Fermi and recoil en- 
ergies, Ek F and e q , respectively, are defined below. This 
density describes both the one-dimensional interacting 
bosons and the noninteracting (degenerate) Fermi gas. 
Since the density becomes spatially modulated after the 
interaction, regions of high density, or focal spots, can 
arise when the spatial harmonics at a certain time su- 
perpose in phase at certain points in space and out of 
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phase in others. In addition, the spatially modulated 
density exhibits a Talbot effect, resulting from the peri- 
odic self-imaging or wave packet revival of each fermion's 
time-dependent, single particle wave function. Simulta- 
neously, the inhomogeneous dcphasing due to the initial 
Fermi distribution of momenta causes a decay in time 
of the spatial modulation. Since the Talbot effect is a 
coherent, single particle phenomenon and the distribu- 
tion of initial momenta represents the interactions in the 
system, in some sense the decay of the modulated den- 
sity in the condensate measures the crossover from single 
particle to many-body physics in an interacting system. 



II. THE MODEL 



We consider N bosons with hard core interactions in a 
one-dimensional "box" of length L with periodic bound- 
ary conditions. The Hamiltonian H for this system cor- 
responds to particles with mass M and radius a — ► 
interacting via effective point or ^-function interactions 
in the limit of infinite interaction strength: 

H = -^Y,Q^2 +2E f ( I *" I j)' With7 ? ^00. 

i=l 4 i=ij 

(1) 

The effect of the interactions is to introduce nodes in the 
many-particle eigenstate wave functions whenever two 
coordinates coincide: the bosons are impenetrable. The 
problem of impenetrable bosons in one dimension has a 
very close relationship to the problem of free fermions. 
The energy eigenfunctions differ only by a multiplica- 
tive function with value ±1 which symmetrizes the anti- 
symmetric fermion cigenstatcs to form boson states, and 
the corresponding eigenvalues are identical @,|. The 
eigenstates can be labeled by a set of wave vectors {ki} 
with i — 1,---,N, ki = -r-Tii, and m integers satisfy- 
ing rii ^ rij. We take N odd. The N even case can be 
treated similarly ||. 

For free fermions the ground state wave function is the 
following Slater determinant of plane waves, 

¥ F) ( Xl X 2 ,---X N )= I det e 2nin jXm /L_ 

VN\L N (l-N)/2< nj <(N-l)/2,l<m<N 

(2) 

The fermionic wave function corresponds to occupying N 
one-particle states with the lowest energy satisfying the 
Pauli exclusion principle. The result is a "Fermi sea" of 
occupied states in the interval k € [—hp, hp] for 

kp/ir =(N- l)/L = Po = 1/d, (3) 

where po is the mean density and d is the interparticle 
distance. The excited states are constructed by emptying 
some states of the Fermi sea and occupying states with 
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|fc| > hp. The corresponding bosonic wave function can 
be written simply as 



det e 2irinji m /I 
(l-7V)/2<rij<(Af-l)/24<m<Af 



(4a) 



1 TT \ e 2 ^n/L _ e 2irzx m /L / 4b) 

VJV.lv i< n<m <N 



Since the bosonic wave function (for all the states of the 
spectrum) is even under permutation of any pair of co- 
ordinates Xi with Xj , Cill equivalent description of the 
ground state wave function is 

*1 B) (X1 < X 2 < ■ ■ ■ < X N ) = -L= det e 2^n ]Xm /L_ 

y<L N {l-N)/2<n J <(N-l)/2S<m<N 

(5) 

In second quantization, the Hamiltonians for fermions 
and impenetrable (or hard-core) bosons are very similar: 



H F =Y2e k c\c k , (6a) 
k 

H B =^e k a\a k , (6b) 



where 



Cfc 



oh 

ikx , 



— j= f dx e lkx c x , a k = — ^= f dx e 
V L Jo V L Jo 

(7) 

£fc = fc 2 /2 is the kinetic energy with M, h — > 1, and 
the field operators c x and a x destroy a fermion and a 
boson, respectively, at site x. Despite their similarity, 
there is a basic difference in the commutation relation 
obeyed by such operators, {c Xl c\,} — S(x — x') and 

[a X) <4'] = 5{x — x'), supplemented by the hard-core con- 
dition (a x ) 2 = 0, where [, ] is the commutator and {,} is 
the anti-commutator. In second quantization the hard- 
core interaction has been absorbed into the definition of 
the field operators which must be used to construct any 
state. This is a key point. As a result of the condi- 
tion (aj.) 2 = 0, the boson ground state cannot be easily 
constructed from the vacuum using the creation and an- 
nihilation operators, a k and a k . Furthermore, we do not 
know how to calculate the state rotation caused by inter- 
action operators (i.e., a pulse) on arbitrary boson states. 
However, the fermion energy eigenstates are well known, 
and the ability to determine the action of operators on 
any fermion state allows us to perform tractable calcula- 
tions. 

The state vectors written in the (real) Fock space take 
the same form for fermions and bosons. For example, 
take a state consisting of three particles, | X ■ X ■ X ), with 
x < x' < x". The non-trivial difference appears in the 
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evaluation of off-diagonal matrix elements. For the pre- 
vious state, 



t t t 

= C i C X > C x> 

\x,x',x")B = aial,al„\$), (8b) 



|z,z',:r">F = 44'i<'|0>> ( 8a ) 
t t t 



the following off-diagonal matrix elements illustrate the 
difference in question (xq < x), 

(x Q ,x,x"\clc x '\x,x',x") F = -1, (8c) 
(x ,x,x"\ala x '\x,x', x") B = +1- (8d) 

This difference arises from the commutation relations. As 
a direct result, the commutation relations allow for the 
possibility of having more than one interacting boson in 
the k = state as long as two bosons are not at the same 
position, whereas for free fermions the exclusion principle 
limits the occupation probability of the k = state to one 
or zero. 

Information about the condensate in this one- 
dimensional model is contained in the (single particle) 
density matrix, 

p(x,x',t)=(*(t)\ata x ,\*(t)). (9) 

For the ground state, the density matrix, 

p(x,x') = (y \aia x ,\y ), (10) 

was analyzed by Lenard The condensate wave func- 
tion ^> c (x) is defined as the eigenstate of p(x, x') with the 
highest eigenvalue N c : J Q dx' p(x, x') ^ c {x') = N c ^ c (x). 

According to the Penrose and Onsager criterion ||, 
a true Bose-Einstein condensate corresponds to N c ~ 
0(N). For the present one-dimensional case with hard 
core interactions, Lenard showed that N c ~ 0{\/N) ||. 
This implies that in the thermodynamic limit the par- 
ticle density vanishes in the condensate. However, for a 
finite system there can still be a large number of particles 
in the condensate. This vanishing density is connected 
to the fact that there is no Bose-Einstein condensation 
in one-dimension for a free Bose gas at finite tempera- 
ture. Even though we are considering here the limit of 
zero temperature, the quantum fluctuations generated by 
the hard core interactions are enough to destroy the con- 
densate in the limit of infinite particle number. Given 
the above proviso, our model is relevant for the under- 
standing of interactions since an infinite, non-intensive 
number of particles remains in the lowest lying state. 

For fermions the ground state wave function corre- 
sponding to the Fermi sea, [Eq. (0)], can be written 
in second quantized form as the state vector 



\f)= n 4i0>: ( n ) 



and the ground state density matrix, Eq. ( 10) with a x 
c x , can be evaluated easily as 
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PF (x,x') = I J2fke ik( - x '- x) . (12) 
fc 

The Fermi factor fk satisfies fk = 1 for \k\ < kp and zero 
otherwise. From Eq. ( |l2| ) the eigenstates of the fermionic 
density matrix are plane waves e iqx with eigenvalues f q . 



F),is 



The ground state energy, given by {F ^2 k £kc\.Ck 

Ef ] = E { B) ~ N(hk F ) 2 /(6M). 

We are interested in the time dependence of the boson 
state vector \&(t)) and density, p{x,t) = p(x,x,t), after 
a pulse is applied with the form of an impulsive, periodic 
phase grating, 

N 

H'{t) = -2XS(t) cos q Xi . (13) 

i=l 

One can prove [see Appendix [X| that the action of this 
pulse on either the fermion or boson ground state leads 
to the same time-dependent, spatially-modulated den- 
sity at zero temperature given the above connection be- 
tween fermion and boson eigenstates and eigenenergies. 
The crucial feature in the proof is that both the pulse 
operator and the diagonal elements of the density ma- 
trix depend only on the boson coordinates (the position 
x and the diagonal operator a].a x in second quantiza- 
tion). These operators preserve the symmetrization of 
the fermion states to form boson states. Therefore, while 
the occupation probability of A;-states is distinct for the 
interacting bosons and free fermions, the density of the 
Fermi system mirrors the collective behavior of the many- 
body (correlated) boson wave function for particles sep- 
arated on average by a distance d. Furthermore, the fact 
that the fermions act as free particles enables us to cal- 
culate the results for the interacting Bose gas exactly. 

^From now on we will consider the pulse acting on a 
system of fermions, the complete Hamiltonian in second 
quantization being 

H = - A<5(t) (p, + p\) =H - XS(t) ( Pq + p\) , 

(14) 



k 



where 



P. 



t \^ 



c\c k + q and p q = pl q . (15) 



Note that we are omitting the subscript F for fermions 
and that the pulse Hamiltonian in Eq. ( |l4| ) is the second 
quantized version of H'(t). The following commutation 
relations and identity, derived from {c^, cLj = Sk.k' an d 
{cfc, Ck'}k=jLk' = as well as Eqs. (0) and (^5[), will prove 
useful below: 

[ck,P q ] = Ck-q, (16a) 

[Ck,p\] = Ck+q, (16b) 

C x = ^=J2 C k eikX - ( 16c ) 
^ L fc 
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III. RESPONSE TO THE PULSE 



For times t < the wave function is the Fermi sea, 

|*(t < 0)) = \F). (17) 

For t = + we have from the integration of the 
Schrodinger equation (|lj) in the impulse approximation, 

|*(t = 0+)) =e iX ( p " +p D\F), (18) 

and for times t > 

=e- iaot e- iX ( p * +p <)\F). (19) 

We are interested in the density p{x, t) from Eq. (^|) for 
positive times, which is given by 

p(a ? ,t) = (*(t)|4a B |*(t)> (20a) 
= ^<*(t)l4c fc ^(t))e^'-^ (20b) 

k,k' 

= ^^^| e -iA(p 9 +pt) c t Cfe , e a(p 9 +4)|^ e i[(fe'-fc)x+(e (! -6 fc ,)t] ) (20c) 
fc,k' 

where the third line follows from the second using 
e lHot c k e~ tH<)t = c k e~ leht . The problem has been reduced 

to simplifying the operator e~ jA ( p?+p 9)c[,c fc /e jA ( p « +p «) . 
Each term in its power series expansion couples the Fermi 
sea to states within and outside of the Fermi sea. The key 
to this calculation is that we can evaluate this operator 
explicitly, summing the effect of the pulsed interaction to 
all orders in A. 

Rewriting e _<A ( p « +p S) 4c fc 'e a ( p " +p ') as 

e _iA( P5l +pt) c t e a(p ( ,+pt) e _ !i A(p ( ,+4) Cfe;e iA( P5l +pt)^ we de _ 

fine 

Cfc(A) = e - a ( p « +p l)c k e a ( p « +p l), (21) 

where the q-dependence of Cfc(A) is implicit. Taking the 
derivative wit h resp ect t o A and using the commutation 
relations (16a) and ( |l6h| ), 

^) = ie-»M) [ Ck , Pq + p\) ^M) 

= i(c k - q {\)+c k+q (\)). (22) 

The solution to this equation with the condition c k (X — 
0) = c k is 

OO 

Cfe(A)= i s Js(2X)c k - sq , (23) 

s— — oo 

where J s is a Bessel function of integer order s. 

The action of the pulse transforms the single fermion 
annihilation operator into a coherent superposition of sin- 
gle fermion operators. In fact, from Eq. (|lj) we can con- 
firm that Cfc(A)e -l£fet is the time-dependent Heisenberg 
annihilation operator. Substituting 



G 



(F\c{(\)c k ,(\)\F) = J2 i ^~'M 2 >)M^)Sk'- a ' q ,k-« l fk'-.' q 

s.s 



into Eq. (20c), letting k — sq — ► k and k' — s'q — ► k\ 



and summing over k' , we find 

POM) = 7^^* s '"V s (2A)J s ,(2A)e I [( s '-^^ + ^+-- e ^'^ t ]/fc. 



ft; s.s' 



(24) 



^From this equation the calculation shows that, as far 
as the density is concerned, each particle with initial 
momentum hk acts independently, multiply-scattering 
the pulsed interaction to form momentum components 
h(k + sq). The components of the perturbed particle 
then propagate to time t with their free energies, £k+sq- 
The total density at t is formed by the interference of 
each particle's individual momentum components, aver- 
aged over the initial momentum distribution fk . This av- 
eraging process determines whether collective effects ow- 
ing to particle-particle interactions dominate (q <C 2kp) 
or whether the coherent interference of each particle with 
itself dominates (q ^> 2fcp). The dephasing caused by in- 
teractions results in a decay of the density pattern as a 
function of time. 

Transforming 2ir^2 k — ► L J dk, we can perform the 
integral by noting the Fermi function fk restricts the lim- 
its of integration to the interval, —kp < k < kp . There- 
fore, 

fcp 

I^/ fc e*K«*+-«-^«)*] = _L e V( S 2 -s' 2 W2 j dk e ik g (s-s')t (25a) 

k , 

= 1 c ic 2 (s 2 -s' 2 )t/2 MkFq(s' - s)t] ^ 
tt q(s' — s)t ' 



and we can make a change of variable in Eqs. (24) and 



(251), s' — s — j, to give 



p( X)t) = lYiij^-Mm ^Ftii] y j s{2X) j {2X)e -isfit_ 
tt z — ' qit — ' 

3 J s 

(26) 

Using the sum rule [M 

JsC2X)J s+3 (2X)e- 2lsv = r J ' J;j(4\ sin <py jv (27) 

S 

for cp = jq 2 t/2 yields the desired result, 

p{x,t) = - Y Jj (4A sin jqH/2y™* sin ^ kp ^ ( 28 ) 
tt ' ' qjt 

j = -oo 

or 

oo 

p(x, t) = p + 2 Jj(AXsin jhq 2 t/2m) cos(jqx)Ij(t). 

3 = 1 

(29) 
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The total particle density is spatially modulated as a 
function of time, and each spatial harmonic decays in 
time according to the decoherence function, 

_ 1 sm[k F qjt] _ 1 s'm[hk F qjt/M] 
iU ~~ 7f qjt ~ n hqjt/M ' [ ' 

owing to the uncertainty in the initial A:- vector within the 
Fermi sea. We have reinserted the H and M to stress the 
time scales of the solution. We now examine two regimes 
of this exact solution, q ^> 2k F (where collective effects 
are negligible) and q < 2k F (where collective effects de- 
grade the density modulation). 



IV. RESULTS 



A. q > 2k F 



For the q ^> 2k F case when hk F qjt/M <C 1 but 
Hq jt/(2M) is arbitrary, we can ignore the decoherence 
by setting Ij(t) — ► Iq = kp/rr = p , the background 
density. In this limit the total density is a periodic func- 
tion of time with period t f — 2tt / {hq 2 / (2M)) . This is 
the Talbot effect from classical and atom optics, where 
the periodicity of a diffraction grating is imposed on a 
transmitted wave, and tt is the Talbot period |Ic|-|T^1. In 
those cases the transverse wave intensity or density in the 
Fresnel regime of diffraction becomes periodic in space 
and time in exactly the manner of Eq. (|2g|). The ini- 
tial, uniform density evolves into a spatially-modulated 
density and back again. 

If 4A > 1, a third time scale, t/ = 2M/{A\hq 2 ) = 
7W(8ttA), is implicit in the expression for the density, Eq. 
( p9| ) . For t ~ Tf the particle distribution is focused to an 
array of focal spots with density peaks of magnitude > po 
at the positions qxf — 2nm for all integers m |l^,|l3| ■ As 
a result of the Talbot effect, these focuses also reappear 
at later times t = nT F /2± Tf for integers n > 0. The 
focal positions qxf are shown in Table 1. 



qx f 


n odd 


n even 


t = riT F /2 + Tf 


277171 — 7T 


2ivm 


t = utt/2— Tf 


27TTO 


2nm — 7r 



Table 1. Focal positions for 4A > 1, q 3> 2k F 

Hence, the spatially-modulated pulse acts on the many- 
body wave function as a periodic array of lenses. 

In Fig. la we plot the normalized particle density, 
p(x,t)/po, versus \qx \ /tt < 1 on the horizontal axis and 
< </tt < 0.5 on the vertical for A = 2.5 and k F — > 0. 
The highest density spots are white while the lowest are 
black. The density p(x,t)/p at t/r T = 0, 0.5 is one. Fur- 
thermore, we can see that the density in this time period 
is symmetric with respect to t/rr = 0.25. The white spot 
near x = and t/T F ~ Tf/vr — 0.016 is the first focus; 
the white spot near x = and t/rr ~ 1/2 — 17 /tt is the 
second. The density for the time period 0.5 < t/vr < 1 
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(not shown) is identical to Fig. la, shifted by half of the 
spatial period. In Fig. 2a the normalized density ver- 
sus time, p{qx = 2irm,t)/ po, is plotted for A = 2.5 and 
kp — > 0. This is the slice of Fig. la along x = 0. The 
Talbot effect is clear as the amplitude fluctuations repeat 
with period tx- 



B. q < 2k F 

This situation corresponds to the limit for which the 
momentum "kicks" due to the applied pulse arc smaller 
than the momentum uncertainty in the ground state, and 
one expects the Talbot effect to be washed out for higher 
harmonics. From Ij (t) we see that the largest harmonic 
that survives up to the Talbot time is that for which 
kpqjTx < 7r, implying j ~ q/{Akp). Higher harmonics 
die away before the pattern can be "reconstructed" at 
Tx- In quasi-classical language, a particle that absorbs 
momentum jhq takes a time r c ~ d/(hjq/M) to travel 
the inter-particle distance d = n/kp and collide with 
its neighboring particle. If r c is much longer than the 
Talbot time tt, the particle acts as a free particle for 
multiple Talbot times, and the response of that harmonic 
will be insensitive to the particle-particle interactions. 
The condition r c ~ tt gives an estimate of the highest 
harmonic j^ ax surviving up to the Talbot time, 

•T ^ 1 /oi \ 

Akp 

which agrees with the consideration kpqjrx < 7r above 
from the exact result. 

To continue the analogy with classical and atom optics 
including the decoherence factor Ij(t), the density 
of Eq. ( p9| ) is isomorphic to the behavior of the atomic 
density when a beam of neutral atoms with beam diver- 
gence 9 — 2hkp jM a v z -C 1 passes through a standing- 
wave phase grating with periodicity 2ir/q, where M a is 
the atomic mass, v z = z/t is the longitudinal beam speed, 
and IhkpjMa plays the role of the transverse velocity 
spread u x . In that case, the probability distribution of 
transverse atomic velocities is taken as P{v x ) = l/u x for 
v x G [—u x /2,u x /2], just as the occupation of /c-states in 
the Fermi sea is P(k) — l/(2kp) for fc S [— kp,kp] in 
one-dimension. A similar correspondence holds for the 
divergence angle, = 2kp/k Zl of a light beam in the 
paraxial approximation passing through a periodic opti- 
cal phase grating. As a result, Ij(t) corresponds to the 
Dopplcr dephasing caused by a uniform, inhomogeneous 
distribution of particles with initial wave vectors between 
±kp with respect to the momentum kick wave vectors 
jq. The sinusoidal behavior of Ij(t) arises from the uni- 
formity of the Fermi distribution over a finite interval 
of /c-states. Other inhomogeneous distribution functions, 
such as a thermal distribution of particle velocities in the 
atom beam case, may lead to a smooth decay in time. 

Again, if 4 A > 1, we can have the situation where 
T f ^ T c r T, implying the particle distribution comes 
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to its first focus before the spatial modulation washes 
out. Actually, the stricter condition, Tf(2jhkp/M) < 

(4A) 3 ^ 4 An/q, is required. This condition estimates the 
highest harmonic contributing to a focus for which the 
spot size for kp — ► 0, w ~ 47r (4A) 3 ^ 4 /q, is not broad- 
ened by the Fermi distribution of momenta (w). As a 
result, we have 

£ ax ^2 7 r(4A) 1 / 4 ^|-. (32) 

Even if q < 2kp, j^ax can be significantly greater than 
one for 2tt (4A) 1/4 > 2k F /q. 

In Fig. lb we plot p(x,t)/po versus \qx\ /tt < 1 and 
t/vp for A = 2.5 again but now kpjq = 0.1. The pure 
Talbot effect is washed out as the modulation harmonics 
of the density clearly damp in time, and the symmetry 
with respect to t/rp — 0.25 is broken. The first focus is 
still apparent while the second is vague. This is a regime 
where the initial focusing occurs for all significant j (i.e., 
imax ~ 28 and J 2 s(4A) <Cl), In Fig. 2b the normalized 
density, p(qx — 2 / nm 1 t)/po, is plotted for this case. From 
Eq. (E9I) the largest significant harmonic surviving to the 
first Talbot time is j^ ax ~ q/{^kp) ^2 — 3. 

To emphasize the dephasing, Fig. 3 shows a compari- 
son between the density at the first (Fig. 3a), second (Fig. 
3b), and fifth (Fig. 3c) focuses (t ~ Tf, tt/2 — Tf, and 
Tp + Tf, respectively) for A = 2.5 and kp/q — 0,0.1,1.5. 
The focal density's peak is diminished and its width 
broadened by the particle collisions as time progresses. 
While we calculated j£ ax ~ 28 for kpjq = 0.1, j£ ax 
is approximately 4 for kpjq = 1.5. The first focus for 
kpjq = 1.5, denoted by the dashed line, is clearly flat- 
tened and broadened in Fig. 3a while the focal density for 
kpjq = 0.1 is indistinguishable from the kpjq = case. 
By the second focal time in Fig. 3b, the kp/q= 1.5 case 
has returned to nearly uniform density (i.e., the t — > oo 
limit), and the kpjq — 0.1 case, denoted by the dash-dot 
line, has started to damp. In Fig. 3c both cases kp /q^Q 
have nearly uniform densities at t = Tp + Tf, the focus 
after one Talbot time. From these plots it is clear that 
a comparison of the normalized density as a function of 
time can act as a sensitive measure of the time scales for 
which many-body effects become important compared to 
single particle effects. 

V. PERTURBATION LIMITS AND 
FUNDAMENTAL PROCESSES 

By looking at the perturbative limit of these results, 
insight is gained into the fundamental processes taking 
place in this system. In particular, the cases of large 
spatial period (q 2kp) and small spatial period (q 3> 
2kp) pulses can be further distinguished. Taking the 
perturbative limit of Eq. ( ]29j ) to first order in A, only 
the lowest-order harmonic survives. The result is 
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p(x, t) ~ po \l + 4A sin (^fj cos qx ^ff j , (33) 
where Vp = hkp/M is the Fermi velocity. This first- 



order density modulation can be traced to Eq. (20c) 

taking e ±iX ( Pq+p ^ « 1 ± iX (p q + p\). Two time scales 
are present in this limit, the Talbot time tt and the de- 
phasing or collision time r c ~ (qvp) -1 . By examining Eq. 
( |33| ) for small times, two different pictures of the quan- 
tum processes emerge, depending on the size of q/2kp. 
Note that, for t/r>r < 1 or t/r c <C 1, the signal does not 
decay. 

For hq 2 t/(2M) < 1, we can rewrite Eq. (|3|) as 

p(x, t) ~ po ( 1 + 2A— cos qx sin qvFt ) (34a) 

V k F J 

= Po ( 1 + A- — [sing [vpt — x] + sinq\vpt + x]] ) . (34b) 

V k F J 

If this expression is to be valid for times qvpt ~ 1 when 
hq 2 t/(2M) <C 1, the condition, g <C 2fep, is necessary. 
This is the lowest-order quantum scattering limit where 
interactions dominate. The density modulation has a 
small amplitude, 2Xq/kp, and is periodic in time with 
frequency qvp. The sinusoidal time dependence can be 
traced to the uniform Fermi distribution of momenta. 
Thus, in the limit q -c 2k p a small amplitude modu- 
lation can be created on a length scale 2ir/q such that 



L^$>2n/q^$>d = L/N.ln some sense Eq. (34a) looks 
like a classical Doppler modulation of the density as a 
function of time. For q <C 2kp the scattering can be 
seen as an effective process where only the small fraction 
(~ q/kp) of particles near the Fermi surface can scatter 
the pulse due to the exclusion principle, and these parti- 
cles act independently on each side of the Fermi surface. 
The pulse creates a superposition of two amplitudes for 
each of the particles, either with wave vectors ~ kp and 
~ (kp + q) on the right side of the Fermi surface, or with 
wave vectors ~ (— kp) a nd ~ (—kp — q) on the left. The 
density pattern of Eq. ( |34b| ) forms as each particle in- 
terferes with itself with the same relative phase ~ qvpt 
but opposite relative momenta, hq and — hq, for the right 
and left running particles, respectively. At longer times 
the recoil energy of the scattering, H 2 q 2 /(2m), emerges 
in Eq. (|3^), signifying the breakdown of this limit and 
the decay in time of the modulation amplitude. 

In Fig. 4a the density, p(x = 0, i)/po, is plotted for A = 
.05 a nd k p j q = 6. For early times the signal is given by 



Eq. ( |34a| ). At later times Hq 2 t/(2M) ~ 1, a beating with 
the Talbot oscillations appears as a result of the product 
of sinusoidal functions in Eq. ([53|). This beating occurs 
only because the Fermi distribution is piecewise uniform, 
giving rise to the time dependence sm[qvpt]/qvpt. Care- 
ful observation shows that the signal depends only on the 
beat frequencies, \2kp/q± l|/ir = ll/ry, 13/tt- 
For qvpt -C 1, we can rewrite Eq. (|33|) as 
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p(x,t) ~ po (l + 4Asin (i^j^j cos qx^j . (35) 

If this expression is to be valid for times hq 2 t/(2M) ~ 1 
when qvFt <C 1, the condition, q 3> 2kp, is required. 
This limit recovers the Talbot effect and is the lowest- 
order quantum scattering limit where single particle ef- 
fects dominate. Note that the amplitude of the mod- 
ulation does not depend on the factor kp. For times 
hq 2 t/(2M) ~ 1 the relative phase from the initial ki- 
netic energy of the fastest particle in the Fermi distri- 
bution is insignificant compared to the phase from the 
recoil energy gained by each particle from the pulse, 
\h 2 (±k F + q) 2 /{2M) - h 2 k F /(2M)\ ~ h 2 q 2 /(2M). As a 
result, this looks like a case of quantum scattering from 
a state having k ~ 0. 

In Fig. 4b the density, p(x = 0,t)/p , is plotted for 
A = .05 and kpjq = .125. The signal rises at early times 
as predicted by Eq. (p5|). At later times qvpt ~ 1, 
the beat frequencies of Eq. fl33|), \2kp/q ± l\/vr = 
5/(4tt), 3/(4tt), are apparent. In Fig. 4c for A = .05, 
we plot the degenerate case, 2kp = q, where the two ef- 
fects oscillate in phase to produce a modulation which 
goes as sin 2 (hq 2 t/(2M)) /hq 2 t/(2M) at x = 0. 



VI. CONCLUSION 

The exact results showing the focusing and Talbot ef- 
fects in a one-dimensional condensate may have direct 
application to a number of neutral atom experiments em- 
ploying confined Bose-Einstein condensates or beams of 
cold atoms from a Bose-Einstein condensate . Hope- 
fully, some knowledge has been gained into the role of in- 
teractions and/or inhomogeneity in causing the decay of 
coherent properties in these type of systems when probed 
by transient interactions. 

This paper has also shown the power of using tradi- 
tional condensed matter techniques both to solve new 
problems in solid state physics and to map problems from 
atomic or optical physics onto solid state equivalents. 
The exact density, Eq. (^9|), is equally valid for a one- 
dimensional Fermi system at zero temperature and could 
be adapted for fermions at finite temperatures. Further- 
more, the correspondence between fermions and bosons 
can be extended within this hard core interaction model 
to include a confining potential in one-dimension, such 
as a harmonic trap. The many-body ground state wave 
function of the bosons is the absolute value of a Slater 
determinant of the first N single-particle energy eigen- 
states of the confining potential. For the oscillator case 
the eigenstates are separated by Hu>, and the condensate 
ground state has the total energy N 2 fno/2. Whether this 
ground state corresponds to a true Bose condensate ac- 
cording to the Penrose and Onsager criterion || is an 
open question. 

In addition, we have shown that the density formed by 
the bosons and fermions after the pulse takes the same 
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form as the density of N independent particles with a 
uniform, inhomogeneous momentum distribution passing 
through a standing wave phase grating. As a result, us- 
ing the analogy with atom optics, we expect that the 
recoil oscillations and focusing effects in the boson and 
fcrmion systems could be recovered with echo techniques 
by applying a second pulse to rephase the different den- 
sity harmonics Q. 
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APPENDIX A: PROOF OF CORRESPONDENCE 
BETWEEN BOSON AND FERMION 
OBSERVABLES 

The relationship between free fermions and hard core 
bosons in one dimension was established by Girardeau 
and Lenard [§J. Each anti-symmetric fermion energy 

eigenstate 4>e (xi...Xn) is a Slater determinant of N 
single-particle states. Each many-particle boson energy 

eigenstate <jr e (x\...xn) has the same energy E e as its 
corresponding fermion state. The boson eigenfunctions 

<t>¥\ Xl ...x N ) = A{x 1 ...x N )<t>f\xi...x N ) (Al) 

are formed by symmetrizing <fii F ' (xi...xn) with the op- 
erator 

N 

A(x\...xn) = Y[ sign(xj — Xi), (A2a) 

A 2 = 1 (A2b) 

The operator A(xi...xn) is a many-particle function tak- 
ing the values ±1 and preserving the sign of the boson 
state under interchange of coordinates. 

The following is a proof using these relationships that 
the (single-particle) density and other operators which 
are diagonal in coordinate space are the same for certain 
states of the Bose or Fermi system. The boson density 
as a function of time is defined as 

p(x,t) — J ... J dXi-.AXN-l^*^ (xx.-.XN-l,X,t)^ B \xi...XN-l,X,t). 

(A3) 

At zero temperature we assume without loss of 
generality that the initial boson wave function 
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v E , ' b - > (xi...xat_i, x, t = 0) is a state prepared from the 

( B ) 

ground state '(xi...xn—x,x) by some many-particle, 
unitary operator U , 

¥ B \t = 0) = u¥ B) = UA(xi...x N - 1 ,x)V i F) (x 1 ...x N - 1 ,x), 

(A4) 



where the second equality follows from Eq. (Al) for the 

( F) 

fcrmion ground state ^ . 

The time-dependent wave function ^/^(t) can be ex- 
panded as a superposition of the boson energy eigen- 
states, 

¥ B ^ ( Xl ...x N ^,x, t)=J2 C e e- iE ^% B) (x 1 ...x N ^ 1 ,x), 

e 

(A5) 

where the time-independent expansion coefficients 



C, 



dx\...dxN-idx(j)* f S B \xi...XN-i, x)^ b \x±...xn-i, x,t = 0) 

(A6) 



are defined by the boson wave function at t = 0. When 
we insert Eq. ( |A5| ) into Eq. (A3), the density takes the 
form 

p{x,t) = Y J J -J dx 1 ..Ax N ^ l Ct,C e e- i ^- E ^ t l h <l ) *} B \x l ...x N 

e,e' 

(A7) 



Using Eqs. Eq. (|Al|), ( |A2| ) and (|A4| ), we can rewrite Eq. 
© as 

P (x, t) = J2l-f dx 1 ..Ax N ^ 1 C*,C e e- i ^ E '- E ^ t / n ^ F \x 1 ...x^ 

e,e' 

(A8) 

and the expansion coefficients of Eq. ( |A6| ) as 

C e = J ... J dxi..dx N -xdx<j)< F) AUAy {E) . (A9) 

The key point of t he pr oof arises here. In order to write 
the density of Eq. (A8) entirely in terms of the Fermi 
system, Eq. ( |A9| ) can not contain the symmetrization 
function A. This implies that A must commute with U , 



[U 1 A(x 1 ...x N -i 1 x)\ = 0, 
giving the coefficients 



dx x ...dx N -xdx(t>< F >U^ 



(A10) 



(All) 



Condition (A10) is satisfied if U depends only on the co- 
ordinates Xi...Xn-i,x. (In particular, the pulse operator 
from Eq. (O) above, 
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u 



exp[i2A^ j cos{qXj)/K\, 



3 



(A12) 



obeys condition ( Ald| ).) If this is t he c ase, t he t ime- 
dependent density is defined by Eqs. (A8) and (All) for 
both Bose and Fermi systems prepared by the operator 
U. 

In general, when Eq. (All) is inserted into (A8), the 



density can be written in second quantized form as 

(F\U^e lHot/h cic x e~ lHot/h U\F). (A13) 



This is Eq. (|20aj) from the text. While the density oper- 
ator in second qu antiz ed form c x c x is automatical ly di 



agonal, from Eq. (A13) it is clear that if condition (AlC) 
is obeyed, the expectation value of any operator which is 
diagonal in x can be calculated using this boson-fermion 
correspondence. 
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FIG. 1. Density Plots, p(x,t)/p , for A = 2.5, 

—0.5 < qx/2-K < 0.5 (horizontal axis), and < t/rr < 0.5 
(vertical axis), (a) k F /q -> 0, (b) k F /q = 0.1. The Talbot 
effect is evident in (a) but washed out in (b). (Grey scale: 
white is high density and black is low density.) 
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FIG. 2. Density as a function of time from Fig. 1 along the 
planes qx — 2nm: p(0,t)/po. (a) kp/q — > 0, (b) kp/q = 0.1. 
In (b) the focusing effect is clearly suppressed at later times as 
higher harmonics damp out more quickly than lower harmon- 
ics. Harmonics up to j^ax ~ 2 — 3 have significant amplitudes 
near t = tt- 
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FIG. 3. Focal density for A — 2.5 over one period, 
—0.5 < qx/2n < 0.5, for Fermi wave vectors, kp/q — (solid), 

0.1 (- • - •), and 1.5 ( ), at times (a) t ~ r/, (b) rr/2 — 77, 

and (c) Tr+Tf. At the first focus (a) the kp/q = 0.1 density is 
indistinguishable from that at kp/q = while the kp/q = 1.5 
case is somewhat damped. By the second focal time (b) the 
kp/q = 0.1 density has begun to damp while the kp/q — 1.5 
density shows little spatial modulation. After one Talbot time 
the focuses in (c) for both cases with interactions are severely 
damped. 
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FIG. 4. Perturbation Theory. The density along the 
planes qx = 2nm from the exact result [Eq. (p9|)], p(Q,t)/po, 
is plotted versus time for A = 0.05. The lowest-order spa- 
tial harmonic dominates the modulated density in all cases, 
reproducing the perturbative result [Eq. (^)]. The beat- 
ing between the Talbot period tt and the dephasing period 
tt I '(2 fcf I ' q) is evident, (a) k F /q — 6 - the density goes like 
Eq. ( pa) ) for times t/r T < 1 but {2k F / q){t / t t ) ~ 1. (b) 
fci?/g = 0.125 - Talbot oscillations dominate according to Eq. 
@ for times (2k F /q)(t/r T ) < 1 but t/rr ~ 1. (c) Degen- 
erate case, k F /q ~ 0.5 - the Talbot and dephasing oscillation 
are in phase with the single time scale tt- In all cases the sig- 
nal begins to decay as the dephasing oscillations beat against 
the Talbot oscillations. 
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